arXiv/ 



Cyclic Operator Decomposition for Solving the Differential Equations 
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We present an approach how to obtain solutions of arbitrary linear operator equation for unknown 
functions. The particular solution can be represented by the infinite operator series (Cyclic Opera- 
tor Decomposition), which acts the generating function. The method allows us to choose the cyclic 
operators and corresponding generating function selectively, depending on initial problem for ana- 
lytical or numerical study. Our approach includes, as a particular case, the perturbation theory, but 
generally does not require inside any small parameters and unperturbed solutions. We demonstrate 
the applicability of the method to the analysis of several differential equations in mathematical 
physics, namely, classical oscillator, Schrodinger equation, and wave equation in dispersive medium. 

PACS numbers: 02.30.-f, 02.60.-x, 03.65.Db 



I. INTRODUCTION 

Various classical and quantum-mechanical problems in 
theoretical physics lead to the necessity of solving the 
linear operator equations for unknown functions and, in 
particular, the differential equations. Exact non-trivial 
analytical solutions of these equations, which include fi- 
nite combinations of elementary operations and special 
functions, are known only for a number of specific cases. 
However, there arc many actual and important cases for 
which such exact solutions were not still obtained even by 
using severe approximations for the corresponding inter- 
action operators. For the cases when exact solutions are 
unknown, some approximate methods are usually used. 
They can be conventionally divided into two types: (i) 
varieties of perturbation theory and (ii) numerical calcu- 
lations (which generally are also based on perturbation 
theory). In spite of significant usefulness and applicabil- 
ity, these methods are not free from various limitations 
and disadvantages. The perturbation theory approaches 
may lead to divergent series, they need sometimes suit- 
able unperturbed solutions, and, finally, they do not pro- 
vide even estimations for precision in most cases (see [l[ , 
and references therein). On the other hand, numeri- 
cal schemes, which arc from the very beginning approxi- 
mate, usually also do not give reliable estimations for the 
precision (some reasonings can be found in 0). More- 
over, they can hardly give an asymptotic behavior of the 
solutions at infinity. Thus, the development of the gen- 
eral method which allows to overcome some of the above- 
mentioned difficulties is the main object of our study. 

In this manuscript we develop an approach based on 
the theory of Cyclic Operator Decomposition (COD), 
which gives the opportunities to obtain solutions (ex- 
act or approximate) of the differential equations with ar- 
bitrary operators. The particular solution can be rep- 
resented by the infinite cyclic operators series, which 



acts the previously determined generating function. The 
cyclic operators and the corresponding generating func- 
tion (COD components) can be specified through the 
given operators in the differential equation. Under the 
convergence requirement, these COD components can be 
chosen in different ways depending on the certain prob- 
lem statement. The procedure differs from the using of 
Born series (or corresponding Neumann series) in the 
perturbation theory [^-|6[ and S-matrix theory of Heisen- 
berg, Feynmann and Dyson [7j . It can be understood eas- 
ily by studying, for example, the difference between the 
formal definition of the generating function and Green's 
function (the last one is derived in some cases by using 
the operator resolvent formalism) [H, 0, [jl-tioj . Generally, 
the proposed series does not require any small parameters 
or unperturbed solutions for the convergence. But, as a 
matter of fact, the procedure can be transformed, under 
some certain choice of COD components, to the "stan- 
dard" perturbation theory with small parameters. For 
the potentials without strong singularities, with reason- 
able choice of the cyclic operators and generating func- 
tion, the corresponding scries usually has uniform conver- 
gence. Some additional features and advantages of our 
approach for analytical and numerical solving the differ- 
ential equations are demonstrated in sections below. 



II. THEORY OF CYCLIC OPERATOR 
DECOMPOSITION 

Let us start from the general case of operator equation 
for unknown function: 



Dip = 0. 



(1) 



Here D is an arbitrary given linear operator and ip is 
an unknown function, which can be a vector or matrix 
of arbitrary dimensionality. This equation can lead in 
particular cases to arbitrary linear differential equations, 
which are considered in examples below. 
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Let us consider a pair of operators G and V, which are 
determined by the following condition: 



D = G—V. 



(2) 



Since the choice of this pair is partly optional, we im- 
pose additional conditions on the operator G: 

3G" 1 , that: GG -1 = I; (3a) 

3ip g j£ 0, that: Gi\) g = 0, (3b) 

where I is the identity operator. Any function ip g , which 
satisfies Eq. (|3b[) , will be called generating function. Now 
we can write the following equation: 



I-G- 1 V)ij = ip 



(4) 



As we can check, under the above-mentioned conditions 
for G and tp g , any solution of Eq.® fulfills Eq.([TJ). Equa- 
tion (Q| can be solved in terms of the following Cyclic 
Operator Decomposition: 



I + &- 1 V + &- 1 VG- 1 V 



(5) 



corresponds to the following particular determination of 
the inverse operator in terms of COD: 



I 



G- 



(8) 



Important feature of the proposed theory is that, while 
the conditions Eqs. (f3a|) and (|3b|) should be fulfilled and 
the convergence of the series is necessary, we still have a 
great freedom of choosing G and corresponding ip g . Gen- 
erally, it gives us opportunities to obtain all the possible 
solutions of Eq. ((XJ) . Sometimes we can naturally choose G 
and ipg in accordance, for example, with the correspond- 
ing initial conditions for Cauchy problem or boundary 
conditions for boundary-value problems. 

In some cases, the exact solution Eq.([5]) can be used 
naturally for obtaining the approximate solution with fi- 
nite number of terms. It can be done, for example, when, 
starting from certain number n, the following conditions 

are satisfied: 



G~ l Vj Aj » [G-'Vj V 



These are the sufficient conditions enabling one to derive 
the approximate solution with the prescribed accuracy. 
Further, the proposed method provides another advan- 
tage if one performs numerical calculations. According 
to the exact solution Eq. (JSJ) , we can use recurrent rela- 
tions when calculating numerically the approximate so- 
lutions. In this case, the calculation of any next term in 
the corresponding series does not require more numerical 
resources than the calculation of the previous one. 



This is the exact particular solution of Eq.(fT]) with cor- 
responding particular COD components determined by 
Eqs.Q and ([3]). The solution makes sense only if the 
obtained series is convergent. This can be achieved in 
different cases depending on G~ X V and ip g , for example, 
if we work in Banach space and corresponding operator 

norm is G~ X V < 1. The convergence of some similar 
operator series was considered also in [4|. 

The theory can be easily generalized also to the case 
of the equations with given sources: 



Dip = ip, 



(0) 



where an arbitrary given function <p describes the arbi- 
trary sources. The unknown function ip could be found 
naturally if the inverse operator is known: ijj = D~ X Lp. 
However, the inverse operator D^ 1 can not be easily 
found for a number of problems. Then, we can write 
a solution of this equation analogously by using COD: 



III. EXAMPLES 

In this section we apply the proposed theory of Cyclic 
Operator Decomposition for the various cases of differen- 
tial equations. Let us first consider the Cauchy problem 
for the equation of classical oscillator. Note that this 
equation, if written in other variables, is the stationary 
one-dimensional Schrodinger equation with given energy, 
and it can be transformed also to the Riccati equation 
by using logarithmic substitution. 

/ + W 2 (i)/-0, 

/(*-) - a, (9) 
f(t b ) = b, 

where f(t) is an unknown function, uj 2 (t) is an arbi- 
trary time-dependent frequency and a, b are arbitrary 
constants. Here, it is natural to choose the components 
for COD as follows: 



i'g+G V 



(7) 



In contrast to the case of Eq.([T]), now we can choose 
ijj g = for some non-trivial particular solutions. Then 
we can obtain the particular solution of Eq.([B]), which 



G 



G' 



dt 2 ' 



0s 



. + b(t-t a ), 



dn dr 2 , V = -u/(t). 



(10) 



If we fix (by our local convention, which we will use 
below) that we write for brevity the same variable upper 
limit of integration as the integration variable and deter- 
mine the successive integration (step by step from right 
to left), we can write a simple expression for the solution: 
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f(t) = [a + b(t - t a )] - dt uj 2 (t.) [a + b(t - t. a )} dt + dt uj 2 (t)dt dt u 2 (t) [a + b(t - t a )} dt 



(11) 



It is important to note now, that the presented series 
(Eg. pip) has rapid uniform convergence at least in any 
interval, where ui 2 (t) is bounded. For example, in the 
limited interval [0, t] the rate of convergence for f(t) can 
be estimated in the following way (we assume here for 
simplicity, that t a = tb = 0, b = 0, and the maximum of 



' (t) in the corresponding interval is C„ 



(n-th term of scries) < \a\ 



(Crnax) t 

(2n)l 



(12) 



In the same way, the convergence can be demonstrated 
for other different COD's, when the cyclic operators are 
bounded for the given generating functions in the given 
relevant interval. Moreover, if additionally Lo{t) is a real 
function, u> 2 (t) > w 2 (0) > and t > 0, we have a de- 
creasing alternating series for the above example, and we 
can estimate the precision of partial sum of the series by 
the value of the last term. 

Now we focus on some particular cases of u> 2 (t). To 
demonstrate that it is possible to obtain a solution with 
any prescribed precision, we consider a case when 
to 2 {t) = [1 - ±sint] > and t a = t b = 0, a = 1, b = 0. 
Calculation of the first two terms in Eq. fTTj) gives 



/(*) = 



1 - i(< 2 -t + sini) 



+ 5. 



(13) 



By calculating the third term in Eq. lfTTj) . we obtain S < 
0.0273 if we consider t in the interval [0, 1]. 

Sometimes we can find also the asymptotic behavior 
of the solution. As an example, we consider the case 
oj 2 {t) = —t a and t a = tb = 0, a = 1, b = 0, where a is 
an arbitrary constant, a > — 1. Using again Eq. (|ll[) we 
obtain exact solution in the following form: 



/(*) = 1 + 



t 



a+2 



(o + l)(a + 2) 



(14) 



(a + l)(a + 2)(2a + 3)(2a + 4) 



By analyzing the corresponding series we can obtain a 
simple upper estimate for the solution (Vt > 0): 



f(t) < 1 + 



a+2 



(o + l)(a + 2) 



exp 



2t2 



a + 2 



(15) 



which gives us the following asymptotic behavior at 



f(t) oc exp 



2n 



a+l 



(16) 



In this case, the same asymptotic can be found also from 
WKB theory (see, for example, quasiclassical approxima- 
tion in (HO). 

Let us now demonstrate the selective choice of 
cyclic operators. For that we consider stationary one- 
dimensional Schrodinger equation (we use below the 
units where H = 1 , m p = 1) : 



${x) = 0, 



(17) 



where ip{ x ) is an unknown function, which describes 
quantum state with energy E in the continuum; A, (3 are 
arbitrary real constants. Without loss of generality (one 
can use scale transformations of Eq. (|17[) ) we can assume 
2E = ra 2 ,f3 = 1. To find the solution of this equation, 
we can choose the components for COD in different ways. 
For example, if one interests in the behavior of ip(x) near 

.2 

x = and in small values of E, he can choose G = 
and use nearly the same technique as in Eqs.(|10p and 
(jlll) . However, this choice can be inconvenient for the 
analysis of the long-range behavior. Another variant of 
choosing the components for COD is the following (we 
use also Eq.© for the particular determination of G _1 ): 



G = rrf 



dx 2 ' 



i> g {x) =C 1 e lmx + C 2 e- 
G- 1 



k=l 



G j 



1 



X X 

dx J dx, 

oo —oo 
2 



(18) 



V = -m 
V = A-e x , 

where C\ and C2 are arbitrary constants. To obtain the 
general solution, we consider the particular case of gen- 
erating function ip g (x) = e mx and corresponding partic- 
ular solution ijjp(x). Then we derive by using a rule of 
infinite geometric series: 



i+£(<Vvc 



k=l 



(1 + im) 2 



1 + S 71 



k=l 



(1 + im)' 



(1 + im) 1 



(1 + 2im) 

(19) 
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From here we obtain 



1 + An e nx ■ ] 



1 



-'■-*■ (k 2 + 2imk) 
n=l fc=l v ; 

and the general solution in the following form: 



(20) 



(21) 



The corresponding series converges at any A and real to. 

In a similar way we can obtain solutions for multi- 
dimensional equations. Let us consider the stationary 
Schrodinger equation with potential surface U(r), multi- 
dimensional Laplace operator A, and energy E: 



[A + 2(E - U (r))] ip(f f ) = 0. 



(22) 



Then we can choose G = A, V = —2(E — U(r)), and ip g 
is any solution of Aif) g = 0. From corresponding COD 
we can obtain a solution: 

ip(r) = [l + A" x y + A^VA^V + ...] Tp g = 0. (23) 

The inverse Laplace operator can be written, for example, 
as A -1 = — F _1 /c _2 F, where F is the Fourier transform 
operator and k is an absolute value of the wave vector in 
this transform. 

Another choice of COD components can be better for 
the finding of the bound states with E < 0. We can 
choose: G = 2E + A, V = 2U(r), ip g is any solution of 
\2E + A]ip g = 0, and write the inverse operator G^ 1 for 
COD as follows: 



G" 1 = {2E + AY 



F~ 



1 



2E-k 2 



F. 



(24) 



In this way, we can calculate in some cases the terms in 
the corresponding COD by evaluating the poles at imag- 
inary values kp = ±i\/~2E. 

Now we consider time-dependent three-dimensional 
Schrodinger equation to demonstrate other applications 
of the proposed method. Let us consider propagation of 
charged particle with arbitrary electromagnetic interac- 
tions (below A(r,t) is the vector potential, and c = 1): 



i^- t -\{iV + A(r,t)) 2 -U(f,t) 



^(f,t)=0. (25) 



Here, we can choose the components for COD in a variety 
of ways depending on peculiar properties of the interac- 
tions. One special choice is the following: 



G 



G- 



V 



d 



ipg = ip(f,t ) = *p (r), 



It 



-i J dt, 

to 



(26) 



iV + A(r,t)) +U(f,t), 



where V corresponds to the time-dependent Hamilto- 
nian. It gives the following solution: 



ip(f,t) = 1 + H) / dtV 



(27) 



+ {-iy J dtV J dtV + 

to t 



This solution can be useful for the numerical calcula- 
tions, namely, for the finding the propagator P, which 
gives: ip{r,t + At) = Ptp(f,t) + 0(At n ). If we use ad- 
ditionally internal time-ordering, we can transform this 
expression to Dyson series (see 0, HH ) ■ 

Finally, we consider the wave equation for electromag- 
netic waves in dispersive medium. We assume that an 
arbitrarily given operator e(r, t) which describes electric 
dispersion does not depend nonlinearly on the field, i.e. 
we still have linear problem. In this case the equation for 
unknown vector potential A(r, i) is the following (see for 
example [H[): 



with the initial conditions 



A(r,t) = 0, 



(28) 



Ht=o) 



S(r), 




R(f). (29) 



(i=0) 



We can choose the following components for COD: 



iJj a = S(r)+ / dtr\r,t)R{r), 



(30) 

to 

t t 
G- 1 = J dt r^t) J dt, V = V 2 

to to 

Then, we can find the solution, which follows from the 
corresponding COD series: 



A{r,t) 



1 + £ ( / dt *) / dt y2 ) 

n ~ 1 Vo t / 



S(r) + j dtr L (r,t)R(f) 

to 



(31) 



For this solution the initial magnetic field is equal to 
[V x S(r)] and the initial electric field is equal to R(r). 



5 



IV. CONCLUSIONS 



In summary, we propose the theory of Cyclic Opera- 
tor Decomposition, which allows one to obtain particular 
solutions of linear operator equations for unknown func- 
tions. In most cases it is possible to obtain all the possible 
solutions, which satisfy the given conditions. We demon- 
strate by some reasonings and particular examples that 
our approach has the following remarkable properties: 
(1) there is a freedom in choosing the COD components 
depending on the certain problem; (2) there is a rapid 
uniform convergence for most of the considered cases; 



(3) it is possible to find the asymptotic behavior of the 
solutions; (4) in many cases when one is analyzing the 
approximate solution, it is possible to estimate the accu- 
racy; (5) the proposed approach gives good opportunities 
for efficient implementation of numerical calculations due 
to the recurrent relations that can be used in COD. 
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